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spectra of various hadrons and their inverse slopes T* at niT—m = 1.5 — 2 GeV 
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in p+p reactions at ^fs = 200 GeV. Due to (multiple) minijet production T* 



in general increases as a function of the hadron mass. Moreover, the T*(rn) 



systematics has a "discontinuity" at the charm threshold, i.e. the inverse slope 



Q_|' of D-mesons is much higher than that of non-charmed hadrons and even of 

CD : 

the heavier Ac baryon. The experimental observation of this characteristic 
behaviour in Au+Au collisions would indicate the absence of c-quark rescatter- 
ing. In contrast, the assumption of thermalized partons and hydrodynamical 
evolution would lead to a smoothly increasing T*(m), without discontinuity 
at the charm threshold. The degree of collective transverse flow, indicated 
by the slope of the T*(m) systematics, depends strongly on whether kinetic 
equilibrium is maintained for some time after hadronization or not. 
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Experimental data on single-inclusive hadron production in p+p reactions at yfs = 23 — 
63 GeV show nearly exponential transverse mass spectra at low transverse momentum [T|]. 
Moreover, the inverse slopes ("apparent temperatures") are practically the same for pions, 
kaons, protons and their antiparticles, i.e. they do not depend on the hadron mass. The 
observed deviation from this behaviour in nucleus-nucleus collisions has been interpreted as 
a signature for collective transverse flow ||. 

In this letter, we discuss how the inverse slope systematics extends to higher energies, i.e. 
p+p reactions at y/s = 200 GeV, and transverse momenta on the order of a few GeV. In this 
kinematic domain, minijet production and fragmentation gives an important contribution. 

We shall also discuss Au+Au collisions, where minijets might rescatter substantially. 
This, in turn, should reflect in a characteristic mass dependence of the inverse slopes. We 
will confront the predictions of two extreme scenarios: superposition of minijet production 
(and fragmentation) without final-state interactions (as in p+p reactions) versus local ther- 
malization of parton matter undergoing hydrodynamical expansion. 

In p+p reactions at ^/s = 200 GeV, we expect that the hadrons with transverse masses 
on the order of a few GeV are dominantly produced via fragmentation of minijets. As a 
first step, we estimate the transverse momentum distribution of c-quarks at midrapidity 
employing the well-known expression for inclusive single-jet production within perturbative 
QCD (pQCD) in leading-logarithm approximation (LLA) ||, 

d 3 a f / \ / \ s da"^ cZ / \ 

^~d?p ^ PP ~ * ° 5 + ^ = J d x ad%b G (x a , /j 2 ) G (x b , /j 2 ) - — — — 5 (s +t + u— 2m 2 c j . (1) 

s, i, u are the usual Mandelstam variables of the parton-parton scattering subprocess, and 
da a /dt denotes its differential cross section in lowest order of perturbative QCD. For 
simplicity, we take into account only the contribution from the gg — > qq process. This is 
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sufficient to illustrate our point. 

The expression for da / dt that accounts for the finite quark- mass is rather lengthy and can 
be found in the literature, cf. e.g. ||. We therefore omit it here. We assume m c = 1.5 GeV, 
Aqcd = 300 MeV, and evaluate the strong coupling constant at the momentum transfer 
scale given by the transverse mass of the produced quark, Q 2 = = m 2 + p\ 0. 

G(x,/j, 2 ) denotes the LO gluon distribution function in the proton, which we take from 
ref. H|. Since we work only in LLA, we employ the same momentum scale in the parton 
distribution functions as in the strong coupling constant, i.e. /i 2 = Q 2 . 

In principle, the hadron spectra could be calculated by convoluting the expression for jet 
production with fragmentation functions ||. However, for transverse masses of a few GeV 
such an analysis would at best be qualitative since, e.g., multi jet production and initial 
state radiation are not included. Also, in this domain the fragmentation functions suffer 
from logarithmic infrared divergences which have to be regulated by a model for soft particle 
production. 

Therefore, to calculate the hadron transverse mass spectra we rather employ the PYTHIA 
event generator 0, using the default parameter settings (version 6.115). PYTHIA simulates 
high energy hadronic and leptonic interactions by implementing a large number of hard 
and soft (sub-)processes, and in particular, a scheme for the nonperturbative hadronization 
mechanism. The model goes significantly beyond pQCD in LLA. It describes not only 
single-inclusive minijet production but also includes multi-jet production and initial state 
radiation. PYTHIA is designed to model the complete event structure (like jet profiles, 
multiplicity fluctuations, various types of correlations etc.) and it has been shown to agree 
reasonably well with experimental observations at collider energies ||. 

The Lund string scheme |J is an integral part of the model used to describe the fragmen- 
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tation of (mini-)jets. They are modeled as one dimensional color flux tubes which decay into 
hadrons: quark- ant iquark pairs tunnel in the color field and the field energy is transformed 
into the sum of the transverse masses mjQ. 

The tunnel probability is proportional to exp(— 7rmy//c), where k is the string tension. 
Thus, the creation of quarks with high transverse momentum is heavily suppressed. As 
a consequence, also the produced hadrons cannot acquire large values of mr- The above 
formula also leads to a strong suppression of heavy quarks. The probability for producing 
a light quark as compared to a charm quark is about 1 : 10~ n ||. The energy E and the 
longitudinal momentum p z of the produced hadrons are determined by an iterative scheme: 
for each hadron the fragmentation function f(z) determines the probability that the hadron 
picks a fraction z out of the available E + p z . The default fragmentation function used in 
PYTHIA reads f{z) ~ ^(1 - z) - 3 exp(-0.58GeV- 2 m^/z). 

The Lund string model has been shown to succesfully describe the nonperturbative 
hadronization in e + e~ annihilation events ||. Moreover, the concept of a color flux tube, 
fragmenting according to a universal fragmentation scheme, has been carried over to hadron- 
hadron interactions. The microscopic models FRITIOF, RQMD, and UrQMD ]10| utilize 
string fragmentation routines for the simulation of soft particle production in p+p, p+A and 
A+A reactions. The only difference to strings from e + e~ annihilation are the leading valence 
(di-)quarks — the remnants of the incident hadrons — as string end-points. The excitation 
of the strings is due to single or double diffractive interactions which can be understood 
and parametrized in the framework of Regge theory (see e.g. ||1TJ). These nonperturbative 



Note that 'transverse' is defined with respect to the string axis. Nothing is said here about the 
orientation of the string with respect to the beam axis. 
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processes account for the major part of the total cross sections at CERN-SPS energies or 
higher, y/s > 20 GeV. Strings which are excited in these processes are preferentially oriented 
along the beam axis of the incident hadrons, leading to much higher typical longitudinal 
than transverse momenta of produced hadrons. 




nijrm (GeV) 

FIG. 1. Transverse mass spectra (at midrapidity, y = 0) of various hadrons in p+p reactions 
at y/s = 200 GeV, as calculated with PYTHIA 6.115. The spectra of the individual hadron 
species include all isospin projections and charge conjugated states. The c-quark spectrum (without 
c-quarks) is calculated within pQCD in LLA. 

In the kinematic region where perturbative minijet production becomes important, the 
color flux tubes may no longer be oriented longitudinally, but according to the pQCD sub- 
process that produces the minijets acquire significant transverse momentum. In the frame 
were the z-axis is parallel to the flux tube, the hadrons are still produced according to the 
above fragmentation function. However, the string axis and the particles' momenta are now 
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rotated with respect to the lab frame. 

Figure [I] depicts the resulting transverse mass spectra. One observes that the PYTHIA- 
spectra follow nonexponential distributions, remnant of the perturbative QCD processes that 
describe the minijet production. Also, the "stiffness" of the spectra increases with the mass 
of the hadron. This will be discussed in more detail below. In particular, the slope of the 
.D-meson spectrum equals that of the c-quarks, which are produced purely by perturbative 
parton-parton scattering^. The reason is that D-mesons can only be produced as the leading 
hadron from a c/c quark jet since the tunneling probability of a c — c pair in a color flux tube 
is practically zero (as discussed above). The m^-distribution of Ac-baryons, on the other 
hand, is slightly "softer" since it involves tunneling of a diquark-antidiquark pair out of the 
vacuum (besides the perturbative production of a c-quark) . 

In Fig. |2], we show the inverse slopes as a function of hadron mass. We compute the 
inverse slope by a fit of the transverse mass spectrum to a Boltzmann distribution, 

1 d 2 N 

— — oc exp(-m T /T*) , (2) 

arriTdy 

We restrict the fit to the range rnx — m G [1.5, 2] GeV. 

As already mentioned in the introduction, in p+p reactions at lower energies, the apparent 
temperatures at small pt were found to be about the same for pions, kaons, protons and 
their antiparticles This changes at higher p T and -Js due to the contribution from 

minijet s. 



2 To obtain the number distribution of c-quarks we have simply divided the differential cross section, 
eq. ([[]), by 40 mb. If we multiplied by two (to include also c-quarks) the quark and D-meson spectra 
would coincide. 
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FIG. 2. Inverse slopes (at midrapidity, y = 0) as a function of hadron mass; PYTHIA 6.115 
predictions for p+p at = 200 GeV, and results from hydrodynamics of p+p and Au+Au 
(calculated on the boundary between mixed and hadronic phase and on the T = 130 MeV isotherm, 
respectively). 

The single inclusive cross section at midrapidity and my — m > 1 GeV is dominated by 
color flux tubes that are no longer oriented longitudinally, but have significant transverse 
momentum (due to the pQCD subprocess). This leads to much higher than in case of 
longitudinally oriented strings. Moreover, the inverse slopes of the rriT—m distributions in the 
lab frame are strongly m-dependent. Such an increase of the inverse slope with particle mass 



can also be extracted from p + p collider experiments [E| at yfs = 540 GeV, and from the 



parametrization of the ^-distributions given in ref . |T3[ . Although that parametrization was 
restricted to the energy region ^/s < 63 GeV, it yields slopes for the pi-, phi-, and D-mesons 
that agree with those obtained from PYTHIA to within 20%, if one simply extrapolates it 
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to v/i = 200 GeV. 

The qualitative reason for this behaviour can be traced back to e + e~ annihilation pro- 
cesses which can be modeled as the fragmentation of a string with fixed energy. The observed 
kinetic energy distribution of produced hadrons show a considerably harder spectrum for 
kaons and protons than for pions [14[ • 

Thinking of the initial state in A+A collisions as a superposition of p+p reactions one 
could attribute any change of the particle slopes to final state interactions. The multiple 
soft pr-kicks that a projectile nucleon can experience while propagating through the target, 
i.e. the Cronin effect, has been shown to be small for pions produced in Au+Au at y/s = 



200 A GeV ||15|| . However, this does not automatically hold true for heavier hadrons. Since 
we have not made any attempt to include such multiple scattering effects, nor modifications 
of the parton distribution functions in nuclei, we restrict the application of the minijet /string 
fragmentation model to p+p reactions. 

In ultrarelativistic head-on collisions of heavy nuclei (A ~ 200), the situation might, 



however, be very different as compared to the p+p case. Transport calculations |16| suggest 
that the initially directed momenta of the partons could be quickly redistributed through 
rescattering. This would then lead to the formation of a very hot (T « 300 GeV) vacuum 
of "macroscopic" size (volume ~ 100 fm 3 ), offering the opportunity to study hot QCD. 

Indeed, if minijets with transverse masses of a few GeV thermalize quickly, energy den- 
sities in excess of 10 GeV/fm 3 can be reached |16| (for Au+Au at y/s = 200 A GeV). This 
energy density is much higher than in p+p since the contribution due to minijets increases 
as v4 2/3 (TJ. 

For definitness we consider Au+Au collisions at y/s = 200A GeV that will be studied at 
the Relativistic Heavy Ion Collider (RHIC) in Brookhaven in the near future. As the extreme 
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case we assume that the produced minijets rescatter so frequently that they thermalize 
locally. The subsequent evolution is then described within hydrodynamics. 

The pQCD processes that produce the minijets in the central region occur on a time 
scale of 1/rriT — 0.1 fm. We assume that one to two rescatterings per particle are necessary 
for local thermalization, and take r = 0.6 fm. This seems also reasonable in view of the 
fact that with To = 1 fm one is able to reproduce the measured single particle spectra for 



central Pb+Pb reactions at y/s = 18A GeV, cf. e.g. |L8 |. At the higher center of mass 
energy of RHIC, the parton density in the central region increases, and therefore a smaller 
thermalization time is expected, cf. also |16| , |l9| . 

As initial conditions we assume a net baryon rapidity density of dNs/dy = 25, and an 
energy density of eo = 17 GeV/fm 3 . Employing the formula of Bjorken ||20| , eqTqtiRt = 
dET(To)/dy, with a nuclear radius of Rt = 6 fm, we obtain an initial transverse energy at 
midrapidity of dE^/dy = 1.2 TeV. We have extracted this value for dE^jdy for central (b < 
2 fm) Au+Au collisions at \fs = 200 A GeV from the minijet /string fragmentation model 
FRITIOF 7.02 Q. It is also compatible with the prediction of HIJING [||]. Note that for 
an isentropic hydrodynamical expansion dEj-fdy decreases with time. On the hadronization 
hypersurface, we obtain dEx/dy = 640 GeV fi~8f . 

The initial (net) baryon density at midrapidity, po, is given by a similar expression 
as eo above, except that dEx/dy is replaced by dNs/dy. These densities are (initially) 
assumed to be distributed in the transverse plane according to a so-called "wounded nucleon" 



distribution, e(r<) = e f(r T ), p(r<) = p f(r T ), with f(r T ) = |y 1 - r\fR\. 

In order to respect boost-invariance, we require the longitudinal flow to have a "scaling 
flow" profile, v z = z/t []2T|2"B|. Cylindrically symmetric transverse expansion ||r8| , |2~I| . |2!) 



is 



superimposed. For T > Tq = 160 MeV we employ the well-known MIT bagmodel equation 



9 



of state, p = (e — 4B)/3, where p denotes the pressure and B the energy density of the QGP 
at T = and vanishing net baryon charge. For simplicity we assume an ideal gas of quarks, 
antiquarks (with masses m u = = 0, m s = 150 MeV), and gluons. 

For T < Tc we assume an ideal hadron gas that includes the complete hadronic spectrum 
up to a mass of 2 GeV. At T = Tc we require that both pressures are equal, which fixes the 
bag constant to B = 380 MeV/fm 3 . The normalization is such that for T — > the pressure 
of the nonperturbative vacuum (i.e. that of the hadronic phase) vanishes. By construction 
the EoS exhibits a first-order phase transition. This "softening" of the EoS in the transition 



region strongly reduces the tendency of matter to expand on account of its pressure [p5|,p6 



For a more detailed discussion of the initial conditions and expansion dynamics in Au+Au 



at RHIC energy please refer to ref. [18 



For comparison, we have also extracted the inverse slopes from hydrodynamics in p+p 
at y/s = 200 GeV. In this case, we employ Rt = 1.18 fm, f{r T ) = Q(R T — r T ), dN B /dy = 0, 
dEx/dy = 2.8 GeV (as obtained from PYTHIA), and (for simplicity) the same tq as for 
Au+Au. 

In Fig. |2| we compare the PYTHIA predictions for p+p with those of hydrodynamics of 
p+p and Au+Au. Within the hydrodynamical solution for Au+Au, strong collective flow of 
quark-gluon matter]] Doppler-shifts T* far above the real emission temperature Tc- If kinetic 
equilibrium in the hot hadron gas is maintained for some time after hadronization (say, until 
the temperature drops to 130 MeV), the collective transverse flow can increase even further, 
and T* ~ 400 MeV can be reached for the charmed hadrons (cf. open squares in fig. 0). 



3 The average flow velocity on the phase boundary to purely hadronic matter is approximately one 
third of the velocity of light 18]. 
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Hadrons produced in the expanding QGP can thus reach comparably "stiff" m-p-spectra as 
those produced in p+p via minijets. 

One also observes that in hydrodynamics T* is nearly proportional to m. In particular, 
there is no jump in T* at the charm threshold and the inverse slope of the Ac is larger than 
that of the D, unlike in the minijet /string fragmentation model. In a thermal environment 
(without interactions), the mass is the only hadron-specific quantity that enters its momen- 
tum distribution. Thus, if the perturbatively produced c — c pairs equilibrate with the QGP, 
the inverse slope of the D-mesons is significantly smaller than in p+p reactions at the same 
energy per nucleon (cf. also the discussion of (pt)d in refs. [p^ , p7jl , and [p^j for the effect of 
c-quark energy loss on lepton radiation in ultrarelativistic heavy-ion collisions). 

In p+p reactions, on the other hand, the initial energy density eo = 1.1 GeV/fm 3 is much 
smaller than in Au+Au. In fact, for our choice of initial conditions the initial state is not in 
the pure QGP phase but in the phase coexistence region. Consequently, on the hadronization 
hypersurface there is practically no collective transverse flow, and the inverse slopes of the 
various hadrons are similar, and equal to the real emission temperature T = Tq- If freeze- 
out occurs deeper in the hadronic phase, e.g. on the T = 130 MeV isotherm, a small flow 
is created due to rescattering in the purely hadronic phase. The pressure in p+p reactions 
at -y/i of a few hundred GeV can not exceed that of the nonperturbative vacuum by far, 
and we thus find no significant transverse expansion p9| |. Therefore, in contrast to the 
minijet /string fragmentation model hydrodynamics can not reproduce the experimentally 
observed |TI| bending of the distributions and the increase of T* with m. 

In summary, we have shown that in p+p reactions at high energy the inverse slopes T* 
of the transverse mass spectra (at midrapidity and tut — m of a few GeV) of various hadrons 
are correlated to their mass. This is due to the underlying pQCD subprocess (i.e. minijet 
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production) that produces fragmenting color strings that are not parallel to the beam axis. 
T*(m) shows a very strong increase at the charm threshold, i.e. the inverse slope of the 
D-mesons is much higher than that of the f2-baryons. 

If a dense QGP is created in central Au+Au collisions, in which u-,d-,s-, and c-quarks, 
and the gluons (up to a few GeV of p?) equilibrate kinetically, the inverse slope of the 
D-mesons decreases substantially as compared to the p+p case, and that of the f2-baryons 
increases. Collective transverse flow of such a hypothetical quark-gluon fluid would estab- 
lish a nearly linear relationship between the inverse slopes and the hadron masses. Thus, 
the T*(m) systematics at mr — m ~ 1 — 3 GeV provides an opportunity to experimentally 
determine the degree of heavy quark rescattering and equilibration in relativistic heavy ion 
collisions. 
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